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1 .  Introduction 


Since  Landau  and  Khalatnikov' s  work^  on  transport  coefficients,  i.e., 

the  thermal  conductivity  and  viscosities  of  liquid  helium,  Khalatnikov  and 

2) 

Matveev  reexamined  the  transport  coefficients  due  to  contributions  by 
phonon-roton  scattering  and  noted  that  the  results  of  Landau  and  Khalatnikov' s 
work  with  a  normal  dispersion  deviate  from  experimental  data  by  as  much  as  a 
factor  of  two.  The  correct  anomalous  dispersion  indicates  three-phonon 
processes  (3PP)  to  be  important  near  absolute  zero  temperature,  and  should 
thus  be  included  in  calculations  of  transport  kinetic  coefficients. 

It  is  well  known  that  the  kinetic  equation  with  temperature  gradient  for 
the  thermal  conductivity  vanishes  for  a  pure  phonon  gas,  and  thus  the 
corresponding  thermal  conductivity  becomes  zero.  However,  the  kinetic 
equation  for  the  first  viscosity  does  not  vanish  for  a  pure  phonon  gas  near 
zero  temperature,  and  3PP  become  important.  The  purpose  of  this  paper  is  to 
evaluate  the  contribution  to  the  coefficient  17  ^  of  first  viscosity  from  3PP 
which  were  not  included  in  the  work  of  Landau  and  Khalatnikov.  For  small 
momenta  the  three-dimensional  dispersion  relation  has  been  evaluated 
explicitly  by  adopting  a  soft  potential  to  give3^ 

«<P>  -  CQp(l  +  7p2  -  5p3  +  ...]  ,  (1.1) 

where  Cq  is  the  velocity  of  sound  at  absolute  zero  temperature.  For  bulk 
liquid  helium  of  density  2 . 18  x  10‘2  A"3,  we  can  use  the  following  values4^ 
obtained  to  fit  the  experimental  data3^  of  the  excitation  energy  spectrum: 


In  Sec.  2  we  consider  phonon-phonon  interactions  in  bulk  liquid  helium 
via  second  quantization  to  obtain  the  characteristic  time  for  3PP.  We  then 
solve  the  kinetic  equation  for  the  coefficient  of  first  viscosity  in  Sec.  3. 
Finally,  in  Sec.  4  we  discuss  numerical  results. 


2.  Phonon-phonon  interaction  in  a  ^rTf-<fiTnni ona]  quantum  liquid 

To  consider  the  phonon-phonon  interaction  in  a  three-dimensional  quantum 
liquid,  we  begin  with  the  Hamiltonian  H  and  the  internal  energy  E(p)^  for 
superfluid  liquid  ^He,  given  by 


H  -  J  d3r  [jv*pv  +  E(p)  ]  , 

E(p)  -  E(pQ)  +  J  d3p  ne (p) 


(2.1) 


(2.2) 


where  v  and  p  are  the  velocity  and  density  operators  expressed  in  terms  of  the 
annihilation  operator  a^  and  creation  operator  a^,  n  is  the  number  density  of 
excitations,  and  p  is  given  by  pQ  +  p',  the  constant  and  variable  parts, 
respectively.  Expanding  the  internal  energy  E(p)  as  a  function  of  p'  and  then 
expressing  this  expansion  in  terms  of  annihilation  and  creation  operators,  we 
obtain 


H  -  Hq  +  H3  +  H4  +  .  . 


^  -  Wa  +  , 


(2.3) 
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where  the  other  six  terms  to  be  included  in  Eq.  (2.4)  are  given  by 

combinations  of  a  ,a  ,  a  and  ,  a^  ,  a^  ,  and  fourteen  additional  terms 
P1  p2  p3  Pi  P2  p3 

with  the  combinations  of  four  creation  and  four  annihilation  operators  appear 
in  Eq.  (2.5) . 

The  nonzero  matrix  elements  of  a  and  a^  are 
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The  matrix  element  <n  |Hn|n  >  is  )(n  +t)J<w,  where  w  -  c  p  and  HA  is  thus  the 

p'U'p  [_  p  2  or  0 

P 

Hamiltonian  of  a  simple  harmonic  oscillator.  and  are  the  anharmonic 

perturbative  Hamiltonians.  Transitions  between  different  states  occur  due  to 
^3’  ^4  anc*  other  higher-order  perturbed  Hamiltonians.  The  transition  matrix 
element  from  state  | i>  to  state  | f>  in  second-order  perturbation  theory  is 


given  by 


<f  |  H  |  i>  -  <f  |  H3  |  i>  +  <f|Hji>  +  ^ 


<f|HJjXj|H.|i> 


E.  -  E. 
J  i- 


(2.7) 


According  to  Eqs.  (2.6)  and  (2.4),  the  first  term  in  Eq.  (2.7)  contributes 

only  to  3PP,  and  the  second  and  third  terms  are  related  to  4PP. 

A  finite  lifetime  of  an  elementary  excitation  in  a  quantum  liquid 

originates  from  two  processes.  One  is  the  collision  between  elementary 

excitations,  and  the  other  is  the  spontaneous  decay  of  a  phonon  into  two, 

three  or  more  phonons.  Since  the  collision  probablity  tends  to  zero  at  very 

low  temperatures,  the  collision  process  will  become  less  important  and 

negligible.  At  very  low  temeperatures ,  below  0.6  K,  the  main  excitations  in 
4 

liquid  He  are  phonons.  Therefore,  we  may  consider  the  decay  of  the  phonon 

8 ) 

into  two  phonons  only,  i.e.,  3PP,  consisting  of  two  processes:  the  direct 
process  of  emission  of  a  phonon,  p^,  p^  -  pj  +  p-j,  and  the  reverse  process  of 
absorption  of  a  phonon,  p^,  P2  +  P3  “  differential  decay  probability 

per  unit  time  in  three  dimensions  is  given  by 
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Since  the  matrix  elements  in  both  processes  are  given  by 
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where  n^  is  the  distribution  function  of  phonons  with  momentum  p,  the  total 


decay  probabilities  per  unit  time  of  the  two  processes  are 
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where  u  is  the  Grdneisen  constant  given  by  (po/Cq) (3CQ/3po) . 

The  collision  integral,  which  is  the  difference  between  the  differential 
decay  probability  per  unit  time  for  the  direct  and  reverse  processes,  can  be 


written  as 
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where  5n  -  n  -  n  with  n  as  the  distribution  function  of  phonons  with 

p  p  p  r 

r  ro  *o 


momentum  p  in  the  equilibrium  state.  Since  the  collision  integral  J^pp(n)  is 
given  by  - (n-no)/r3pp,  we  have 
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The  three-phonon  processes  are  small-angle  processes,  and  thus  p^  and  p3 
are  much  larger  than  p^.  Hence  we  may  use  the  following  relations  to  evaluate 
Eq.  (2. 14) : 


3n 


P1 


L)  w 


(2.15) 
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dw, 
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-WCOS0 


(2.16) 


Substitution  of  Eqs.  (2 . 15) - (2 . 16)  into  Eq.  (2.14)  and  integration  of  Eq. 
(2.14)  over  momentum  space  yields 


*3(u+l)2k4T4 
« - — — [1 
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3 .  Viscosity 

In  this  section  we  shall  evaluate  the  coefficient  of  first  viscosity. 
We  first  consider  the  macroscopic  motion  of  liquid  helium,  which  does  not 
depend  on  time,  with  velocity  v  in  the  z -direction  and  the  gradient  directed 
along  x-axis.  Then  the  kinetic  equation  becomes 

v*Vn  -  v  4s  -  J(n)  (3.1) 

x  ox 

The  phonon  distribution  function  in  the  motion  of  liquid  helium  with  velocity 
v  is  given  by 
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[exp( 


Cp-v«p 
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(3.2) 


where  C  is  the  phonon  velocity  (equal  to  Cq  at  zero  temperature). 
Substitution  of  Eq.  (3.2)  into  (3.1)  yields 


n  (n  +1) 
o  o 


_C£  £y 

kgT  ax 


cosflsinflcos^  -  J(n) 


(3.3) 


where  here  6  is  the  angle  between  p  and  the  z-axis.  Since  4PP  establish 
energy  in  the  equilibrium  state  and  the  distribution  function  of  phonons  is 
related  to  the  variation  of  phonons ,  which  characterize  the  chemical  potential 
a,  the  distribution  function  becomes 


i* 


Expanding  Eq.  (3.4)  as  a  function  of  a'  to  obtain  the  deviation  of  n  from  nQ 
and  considering  the  linearity  of  the  left-hand  side  in  Eq.  (3.3),  we  can  write 


5n  -  n  -  nQ  -  -  no(nQ+l)a'  ,  (3.5) 

a'  -  acos0 sintfcos^  (3.6) 

Replacing  J(n)  by  J^pp(n)  in  Eq.  (3.3)  and  integrating  over  momentum  space,  we 
obtain 
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costf sin0cos^(— — ) 
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Since  the  collision  integral  33pp(n) 
use  Eqs .  (3.5)-(3.6)  together  with  Eq. 
of  collision  integral: 


can  be  replaced  by  * (n'n0)/r3PP>  we  can 
(2.17)  and  w  -  VCQ/W  for  t*ie  evaluation 
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From  comparison  of  Eqs.  (3.7)  and  (3.8)  we  have 


a 


30 X4p  C5 
o  o 

x3(u+l)2(kBT) 


_  S2L 
5  3x 


(3.9) 


The  coefficient  of  first  viscosity  is  obtained  from  the  evaluation  of 


the  stress  tensor  associated  with  3PP  given  by 
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Performing  the  integral,  Eq.  (3.10)  becomes 


3!t(3)  (kBX)  flv 
XZ  30*2>(3  Co3  3x 


(3.11) 


We  can  get  the  coefficient  of  first  viscosity  by  comparing  Eq.  (3.11)  with  the 
first  expression  in  Eq.  (3.10): 

3!f(3)KP  C2  _ 

n  h(T) - r - J  -  1.904  x  10"7  J  .  (3.12) 

pn  ir  kgCu+l)^  1  T 

We  note  that  the  coefficient  of  first  viscosity^ ’ 9^  at  temperatures  below  0.9 
K  is  given  by 


*jph(T)  -  3.50  x  10'9  T'1/2  e  /  B  (1  +  2.15  x  10*5T9/2e  ^  )'1  ,  T  <  0.9 


K  . 
(3.13) 


4.  Numerical  Results  and  discussion 

In  previous  sections  we  have  evaluated  the  lifetime  and  the  contribution 


of  first  viscosity.  At  relatively  low  temperatures  (<  0.9  K) ,  Eq.  (3.13) 
reduces  approximately  to 


11 


rj  ,  -  1 . 5 

'ph 


T 


However,  as  Che  temperature  approaches  absolute  zero,  the  roton  density 
becomes  zero,  and  thus  the  contribution  from  the  roton-phonon  scattering  to 
viscosity  will  vanish  so  that  the  main  contribution  comes  from  the  3PP . 
Therefore,  the  temperature  dependence  is  changed  from  T  ^  to  T  Figure  1 
illustrates  the  coefficient  »j  ^(T)  given  by  Eq.  (3.12).  The  coefficient 
increases  slowly,  and  near  zero  temperature  it  increases  very  rapidly  and 
finally  diverges.  In  conclusion,  we  point  out  that  in  consideration  of  the 
3PP  originating  from  the  anomalous  excitation  spectrum  at  low  momenta,  the 
temperature  variation  of  first  viscosity  depends  on  T  \ 
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